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Valleytronics is rapidly emerging as an exciting 
area of basic and applied research. In two dimen¬ 
sional systems, valley polarisation can dramatically 
modify physical properties through electron-electron 
interactions as demonstrated by such phenomena as 
the fractional quantum Hall effect and the metal- 
insulator transition. Here, we address the elec¬ 
trons’ spin alignment in a magnetic field in silicon- 
on-insulator quantum wells under valley polarisa¬ 
tion. In stark contrast to expectations from a non¬ 
interacting model, we show experimentally that less 
magnetic field can be required to fully spin polarise 
a valley-polarised system than a valley-degenerate 
one. Furthermore, we show that these observations 
are quantitatively described by parameter free ab 
initio quantum Monte Carlo simulations. We in¬ 
terpret the results as a manifestation of the greater 
stability of the spin and valley degenerate system 
against ferromagnetic instability and Wigner crys- 
talisation which in turn suggests the existence of a 
new strongly correlated electron liquid at low elec¬ 
tron densities. 


INTRODUCTION 

The valley degree of freedom has a long history 
as a subject of pure and applied research since it 
is an intrinsic property of the band structure of 
silicon and germanium, the historical materials in 
microelectronicJi^. Valley degeneracy had generally 
been viewed as a drawback as it limits the mobility 
of CMOS (Complementary Metal Oxide Semicon¬ 
ductor) devices due to intervalley scattering. Mi¬ 
croelectronics manufacturers have consequently put 
much effort into manipulating valley bands through 
strain, to improve transport properties. This ap¬ 


proach has been successful, and strained silicon has 
been in use in microelectronics since the 90 nm 
nod^. 

More recently, however, the valley degree of free¬ 
dom is becoming recognised as an opportunity, 
rather than a hindrance, and this is leading to the 
emergence of a field of research now known as val¬ 
leytronics in which valleys are exploited in addition 
to charge and spin. Valleytronics has received a re¬ 
cent boost owing to the discovery of graphene and 
other new topical materials also possesing the valley 
degree of freedom and by the proposal of valleytron¬ 
ics device^^Mi^. A vital ingredient to the develop¬ 
ment of valleytronics is valley polarization. Anal¬ 
ogous to spin polarization in spintronics in which 
when achieved under equilibrium conditions, leads 
to key phenomenology such as ferromagnetism, val¬ 
ley polarization can also be expected to yield rich 
and useful physics. 

Experimental research into the physical conse¬ 
quences of valley polarizing a two dimensional elec¬ 
tron system in the steady state is led by studies of 
AlAs and Si based structures. It has been demon¬ 
strated that valley polarisation dramatically affects 
phenomena such as the fractional quantum Hall 
effectfSHni 

and the metal insulator transitior^l^, 
two effects in which electron-electron interactions 
are central. Pioneering experiments performed in 
AlAs indicate that valley polarisation also has a 
strong impact on another effect where electron- 
electron interactions play crucial roles : spin polari¬ 
sation. It has been demonstrated that valley polar¬ 
isation leads to a strong enhancement of spin sus¬ 
ceptibility and symmetric ally, sp in polarisation en¬ 
hances valley susceptibilitj^iSHSi] 

In this article, we firstly confirm the enhancement 
of spin susceptibility by valley polarisation in sili- 


con, which in contrast to AlAs has an isotropic in¬ 
plane effective mass which simplifies interpretation 
of transport phenomena. More importantly, we ex¬ 
plore a new regime in the interaction-disorder pa¬ 
rameter space where a qualitatively new behaviour 
emerges. This is achieved by using electrically con¬ 
trolled valley polarisation in a simple two dimen¬ 
sional electron gas in foun dry co mpatible silicon-on- 
insulator (100) MOSFETiP^ (Metal Oxide Semi¬ 
conductor Field Effect Transistors). We present 
magneto-resistance data which indicate that for low 
enough electron densities, valley polarising the two- 
dimensional electron gas (2DEG) reduces the field 
of full spin polarisation. This represents not only a 
quantitative failure of the single particle picture but 
a qualitative one, in which the observed behaviour 
is opposite to the prediction of the non-interacting 
framework. 




Figure 1. Energy diagram of a 2D electron gas of 
fixed electron density. The diagram depicts the four 
spin-split valley-split subbands. Up (resp. down) arrows 
correspond to spin up (down) electrons and blue (resp. 
red) color corresponds to valley -|- (resp. -). Depending 
on spin and valley polarisation, the system can be either 
(a) valley- and spin-degenerate or (b) valley-degenerate, 
partially spin-polarised or (c) valley-degenerate, spin- 
polarised or (d) valley-polarised, spin-degenerate or (e) 
valley-polarised, partially spin-polarised or (f) valley- 
polarised, spin-polarised. 


RESULTS 

Single particle picture 

Let us first describe briefly how spin polariza¬ 
tion is expected to respond to valley polarisation 
in a non-interacting 2DEG. In the non-interacting 
model, the density of states of a 2DEG is indepen¬ 
dent of energy and can be written as gsg^D^ where 
Dq = mb / 2Trh ?, gs and g^ are spin and valley degen¬ 
eracies, mb is the electron band mass and h is the 
reduced Planck’s constant. In a (001) silicon 2DEG, 
5 s = 5 v = 2 and the system is composed of four in¬ 
dependent spin-valley subbands with equal density 
of states Dq, as depicted in Figure [^a. All states 
are filled up to the Fermi energy Ep = n/gsgvDo = 
u/ADq at zero temperature, where n is the electron 
sheet density. 

Applying a magnetic field parallel to the electron 
gas raises the bottom of the spin down bands com¬ 
pared to that of the spin up bands by the Zeeman 
splitting = gfiBB (Fig. &) where g is the Lande 
g-factor (g = 2 in Si) and /xb the Bohr magne¬ 
ton. Spin down electrons are consequently trans¬ 
ferred to the spin up band: the system spin po¬ 
larises. The spin (valley) polarisation is defined as 
Ps = {n^ — ni)/n (resp. Pv = (n+ — n_)/n) with 
Ti'f and being the spin up and spin down electron 
densities respectively (n+ and n_ are the electron 
densities in the -I- and — valleys). At Pv = 0, full 
spin polarisation is achieved when = 2Ep (Fig.[^ 
c). It follows that the field Bp required for full spin 
polarisation should double when the system is val¬ 
ley polarised because the kinetic energy in the valley 
polarised system (Fig. EJi) is twice that of the unpo¬ 
larised one (Fig. For more details on this single 
particle picture see Supplementary Discussion. 


Experimental determination of Bp 

The field of full spin polarisation Bp of a 2DEG 
can be extracted from measuring the electrical re¬ 
sistance under in-plane magnetic field. Spin polar¬ 
isation has the effect of reducing the ability of the 
2DEG to screen disorder, and as a consequence, in¬ 
creases the resistance due to enhanced scattering 
until, in the simplest case, the magneto-resistance 
(MR) saturates to a constant valuers. This is a sig¬ 
nature that full spin polarisation is reached and the 
spin degree of freedom is completely frozen. 

However, this situation is only rarely observed in 
experiments where it is more common to observe a 
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Fig.j^for a version of Fig. in linear scale). The indicated values for pv correspond to those at full spin polarisation. 
The field of full spin polarisation (plotted as red bullets) is estimated from the field where the MR has reached 97.5% 
of its high-held spin-independent behaviour (grey dotted lines). Values of 5n are given in units of 10^® m“^ in panel 

b. 


shoulder in the magnetoresistanc^i^IllHSlI the 

spin system freezes but the resistance continues to 
change due to spin-independent effects such as the 
coupling of the m agnetic field to the electrons’ or¬ 
bital motioiplJ^. 

In the absence of a comprehensive description 
of the high field behaviour, we follow previous lit¬ 
erature in empirically estimating Bp as the held 
where the MR h as reached 97.5 % of its high held 
dependenc d^^ * ^^ * (The behaviour in the spin polar¬ 
ized regime at high helds was htted to a quadratic 
behaviour with no particular physical meaning. See 
dotted lines in Fig.[^. Results are shown as red dots 
in Fig.[g 

We note that we have also estimated Bp follow¬ 
ing other methods used in the literaturd^ (inhexion 
point in the MR, intersection of the high held asymp¬ 
tote and tangent at the inhexion point etc.). All 
methods provided a qualitatively similar behaviour 
in the changes of Bp with valley polarisation. We 
have chosen the method described above as we be¬ 
lieve it provides a better quantitative estimation of 
Bp compared to other methods which underestimate 
it. 


Evolution of Bp with valley polarisation 

At large density (Fig. and c), we observe an 
increase of Bp with valley polarisation, a behaviour 
which is qualitatively consistent with the single par¬ 
ticle picture (See Refs. and [53] and Methods 
for details on the electrical control and determina¬ 
tion of pv)- It should be stressed, however, that 
quantitatively, the single-particle model fails com¬ 
pletely. The values of Bp are always much lower 
than expected. Also, instead of the doubling of Bp 
from pv = 0 to Pv = 1, we observe only a mod¬ 
erate increase of Bp. For exemple, we measure 
Bp(pv = 0) = 11.75 T and Bp{p^ = 1) = 13.9 
T for n = 2.5 x 10^® m“^ while the single par- 
ticule picture predicts i?p(pv = 0) = 28.85 T and 
Bp{pv = 1) = 51.7 T respectively. These observa¬ 
tions are consistent with the strong enhancement of 
spin susceptibility with valley polarisation seen in 
AlApK^. We note in addition, that the quantita¬ 
tive failure of the single-particle model has also been 
observed in numerous experiments where valley po¬ 
larisation could not be tuned, regardless of valley 

degeneracj*2SH22l 

At lower density (Fig. the single particle pic¬ 
ture fails qualitatively. Here, our data show that Bp 
moves to lower and lower magnetic field as valley 
polarisation is increased. That is, it becomes easier 
to spin polarise a valley-polarised electron gas than 
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a valley degenerate one at low enough density. 
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Figure 3. Exclusion of disorder as the origin of 
the reduction in Bp. a) Bright field scanning TEM 
cross-section of a 16 nm SIMOX device made with the 
same recipe as our sample. The scale bar is 100 nm 
b) High Angle Annular Dark Field scanning TEM im¬ 
age of the quantum well. The image shows the very 
good and similar crystallographic quality of the top and 
buried interfaces. The scale bar is 5 nm c) Electrons 
are pressed against the back interface Sn > 0. d) Elec¬ 
trons are pressed against the front interface for 5n < 0. 
Red dots mark the field of full spin polarisation in both 
panels. Here, n = 1.6 x 10^® m~^ in both panels while 
T = 1.6 K on the panel c and T = 1.7 K in the panel d. 
(This explains the small difference in resistivity at f? = 0 
T and Sn = 0.) The grey dashed lines correspond to the 
high-held behaviour used for the determination of Bp. 


Exclusion of disorder as the cause of the 
reduction of Bp with valley polarisation 

One cause for this observation could be an in¬ 
crease in disorder which is known to reduce 
At first sight, this explanation might seem plausi¬ 
ble since valley polarisation is achieved by pressing 
the 2DEG against the buried Si/Si02 interface (see 
Methods) and hence the 2DEG experiences higher 
disorder due to interface roughness. The ampli¬ 
tude of this effect can be estimated in our device 
where transport can also be investigated at the top 
Si/Si02 interface. At this interface, valley splitting 
is negligibl^^ but disorder is comparable as demon¬ 
strated by our high resolution Transmission Electron 


Microscope (TEM) images (Eig. and independent 
measurements with holes which do not possess the 
valley degree of freedom and show similar mobility 
at both interface^^. This enables us to separate the 
effects of disorder and valley polarisation which are 
mixed at the buried interface. 

Eigures and d show the magnetoresistance of 
the 2DEG for comparable magnitudes of out-of- 
plane electric bias (estimated from the phenomeno¬ 
logical parameter 6n. See Methods) when electrons 
are pressed against the buried or front Si/Si 02 inter¬ 
face. The magnitude of the bare disorder potential 
increases with |i5n|, while valley polarisation is en¬ 
hanced only for Sn > 0.^^ Insight into the variation 
of the bare disorder with dn > 0 is found from the 
comparison of the resistance at (ps = l;Pv = 0) and 
{ps = 0;pv = 1) (highlighted in Eig.|3j:). 

Spin and valley degeneracy can be treated as for¬ 
mally equivalent in the “screening” description of 
transport described in Ref. [231 Therefore, in the ab¬ 
sence of variation of disorder due to the process of 
polarising, one should expect the same increase of 
resistance due to equivalent reduction of screening 
regardless of which degeneracy is liftecP^. 

We find that the resistances are indeed almost the 
same at {ps = l;Pv = 0) and {ps = 0;pv = !)■ 
Therefore, we conclude that the major part of the 
increase in resistance with Sn > 0 seen in Eig. 
at B—0 T can be attributed to a reduced screening 
due to valley polarisation and not to a significant in¬ 
crease of bare disorder. Eurthermore, the increase of 
disorder with Sn > 0 can be estimated to be about 
10% from the observed 10% difference between the 
resistance at {ps = 0;pv = 1) and (ps = 1 ;Tv = 0). 
This estimation of the variation of disorder with |(5n| 
is confirmed by transport at the front interface where 
the entire change in resistance must be attributed 
to a change in disorder. The data at R = 0 T in 
Fig.ii reveal a 14 % increase in the resistance for 
Sn = —0.9 X 10^® m“^ (comparable in amplitude 
to that necessary for full valley polarisation at the 
buried interface). The weak dependence of disorder 
on (5n < 0 is well illustrated in the resistance map 
shown in Eig. where contours of constant resis¬ 
tance run parallel to the constant density lines in 
the relevant regime of density. 

Importantly, Eig. ii shows that the field of full 
spin polarisation is almost unchanged when elec¬ 
trons are pressed against the front interface. That is, 
the 14% change in the bare disorder with Sn in this 
experiments is not sufficient to cause a substantial 
change in Bp. We even note a small initial increase 
of Bp which might be attributed to the removal of 
a small valley splitting at symmetry {Sn = 0) as 
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Quantum Monte Carlo simulations 
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Figure 4. Two dimensional map of the resistance. 

Resistance of the sample as function of front and back 
gate voltages in a log scale. The white line materi¬ 
alises the upper spatial subband edge. The red curve 
represents the pxx = 160014 iso-resistance obtained for 
(n = 10^®m“^,5n = 0 m“^) while the blue one rep¬ 
resents the pxx = 600014 iso-resistance obtained for 
(n = 4.5 X 10^®m“^,^n = 0 m“^). The mobility was 
measured at the white star (see Methods for more de¬ 
tails). 


electrons are moved away from the buried interface. 
Alternatively, this small increase could be due to 
a broadening of the spin band edge due to the in¬ 
creased disorder felt by the 2DEG. Either way, we 
can conclude from this that the increase of the bare 
disorder with |i5n| cannot explain the behaviour seen 
in Fig. leaving valley polarisation as the only cul¬ 
prit. 


As a final confirmation, we have checked that in¬ 
creasing Sn above full valley polarisation does not 
lead to any further change in Bp (see the upper 
curve in Fig. [2^) despite the fact that, as already 
mentioned, the bare disorder potential continues to 
slowly increase with 6n. This has also been con¬ 
firmed in other samples. This set of experiments, 
therefore, allows us to conclude that valley polari¬ 
sation itself is responsible for the reduction of Bp 
observed in Fig. 


The reduction in Bp compared to single par¬ 
ticle expectations is interpreted as resulting from 
electron-electron Coulomb interactions which favour 
spin alignment and this has already been confirmed 

theoreticalljl25h40l 

However, it is particularly chal¬ 
lenging to make quantitatively reliable predictions 
for Bp. Even in the absence of disorder, Hartree 
Fock or random phase approximation (RPA) calcu¬ 
lations do not capture the crucial role of correla¬ 
tions at low density and one has to resort to quan¬ 
tum Monte Carlo (QMC) simulations.^^In our low- 
density intrinsically disordered system, the interplay 
of disorder and interactions further complicates the 
problem. 

Here, we use a Creen’s function quantum Monte 
Carlo approach, known to fully account for the ef¬ 
fect of interactions in the presence of disordeil^lj to 
predict quantitatively the values of Bp in the inter¬ 
acting 2DEC with and without valley polarisation 
(See Methods and ref. HD) . 

The energy per electron for a given spin polarisa¬ 
tion state Ps is given by® 

E = Eo + Eppl ( 1 ) 

We calculate the energy Eq of the spin unpolarised 
and the energy E{pg = 1) of the spin polarised sys¬ 
tem to obtain the energy per electron needed to fully 
spin polarise the system Ep = E{ps = 1) — Eq. The 
predicted field of full spin polarisation is then cal¬ 
culated as Bp = AEp/{gpB)- The resulting Bp as 
function of electron density can be easily compared 
to those extracted from the experiments with no 
adjustable parameter; the only input for the sim¬ 
ulation being the amplitude W of the bare disor¬ 
der (which we can estimate from the peak mobility 
of the sample) and the electron density (See meth¬ 
ods and Ref. for more details on how disorder 
is taken into account in the model). This approach 
has been firmly validated by the successful quantita¬ 
tive comparison to the experimental measurements 
of seven different studies in Si at p^ = 0.®. Here, 
we have extended those calculations to the case of 
Pv = 1 using ppeak = 8000 cm^V“^s“^ measured 
in our sample as a single input parameter for both 
Bp{Pv = 0) and Bp{py = 1) (see Methods for the 
determination of Ppeak)- Assuming that a single in¬ 
put parameter ppeak determined at Pv = 0 is enough 
to describe all the data seems reasonable since we 
have demonstrated that the bare disorder only varies 
weakly with |(5n|. 

The result of the calculation is displayed in Fig- 
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of Bp. (■■••■) Experimental values, (-) non-interacting theory, (—) result of the quantum Monte Carlo simulation. 

b) Ratio of the Zeeman energy at full spin polarisation and Bp as function of rg. Dashed lines correspond to the 
non-interacting case while solid lines correspond to the prediction in the presence of interaction and disorder (The 
calculation are performed for the mobility of our sample), c) Density dependence of ABp = Bp(pv = 1) — Bp(pv = 0). 
Experimental measurements (•) and theory (—). 


ure [^. In both valley-degenerate and valley- 
polarised cases the calculation predicts a much lower 
field of full spin polarisation than the single parti¬ 
cle model. Importantly, the curves corresponding 
to Bp(pv = 0) and Bp(pv = 1) cross so that at 
low density the full spin polarisation is predicted 
to occur at a lower field in a valley polarised sys¬ 
tem. The full dependence of this effect from low to 
high interaction is illustrated in Fig. [Dd. This figure 
shows the predicted ratio of the Zeeman splitting at 
full spin polarisation and the non-interacting Fermi 
energy as function of the interaction parameter 
Tg = l/(7rn)^/^aB (here ub is the Bohr radius) using 
the mobility of our sample. At low rg the behaviour 
is that of a non-interacting system. The effect of 
interaction becomes more and more important as Vg 
increases and the curves corresponding to Pv = 1 
and Pv = 0 eventually cross over at around Vg = 6. 

For a more quantitative comparison between ex¬ 
periments and theory, the experimental values of Bp 
are plotted in Fig. and in Fig. the experimen¬ 
tal difference ABp = Bp(pv = 1) — Bp(pv = 0) is 
plotted together with the prediction. The theory de¬ 
scribes the experiments very well with no adjustable 
parameters, demonstrating that the theory captures 
the essential physics behind the behaviour of Bp. 


DISCUSSION 

As a first remark we should point out that the 
possibility of observing the new behaviour reported 
here results from the cooperative effect of disor¬ 
der and interaction. Indeed, the calculations in the 
disorder free system indicate that electron-electron 
interactions are the leading effect in reducing the 
polarisation energies. Evidence for the crossing of 
Bp{Pv = 0) and Bp(pv = 1) curves is also seen 
for clean systems (See Fig. and Fig. 2 in Ref. 1551 
where it is observed that the spin susceptibility en¬ 
hancement x/Xo{Pv = 1) > 2x/xo(Pv = 0) for large 
enough r^, a feature reminiscent of the crossing of 
Bp curves). However, in Si, this crossing is expected 
to occur at around n = 10^® m“^, too low to be ac¬ 
cessed experimentally in our valley tunable samples. 
By further enhancing the effects of electron-electron 
interactions disorder shifts by a small amount, all 
curves to larger densities, just enough to allow the 
observation of the new phenomenology. The weak 
dependence on disorder implies that our compari¬ 
son is robust against errors in the determination of 
the bare disorder potential. This accounts for why 
Bp is found to be independent of < 0 in Fig. w 
and justifies the use of a single input parameter to 
describe both Pv — 0 and p^ = 1. Nevertheless, in- 


6 














eluding disorder is necessary to achieve quantitative 
comparison. This effect of disorder may explain why, 
in previous experiments in AlAs at similar values of 
Tg, only the enhancement of spin susceptibility and 
not the reduction of Bp with valley polarisation was 
observeeP^. Indeed, in AlAs the mobility was five 
times larger than in our samples and therefore the 
crossing would be shifted to larger r^. However, the 
crossing must have been approached very close since 
the criteria x/Xo{Pv = 1) > ‘^x/XoiPv = 0) is al¬ 
most reached in Fig. 3 of Ref. [121 In addition, we 
should also point out that the comparison between 
QMC calculations and experimental results in valley 
tunable AlAs is complicated by the mass anisotropy 
in the system as concluded in RefHH The situation 
is more simple in narrow AlAs quantum wells where 
effective mass is isotropieP^ and QMC works well.^ 

As a second remark we note that our study fur¬ 
ther confirms that QMC is an appropriate theory to 
predict the polarisation energies of two-dimensional 
electron gasses. Previous studies have shown that 
the disorder free QMC describes correctly the mea¬ 
surements in AlAs (apart from pv = 1 in valley tun¬ 
able AlAs with anisotropic effective mass) and GaAs 
if the finite thickness of the system is includecP^^^. 
We have also shown in Ref. BT] that the QMC in¬ 
cluding disorder describes the available experimen¬ 
tal data in silicon at Pv = 0. Therefore, QMC is 
able to determine the polarisation energies of most 
investigated samples without adjustable parameters. 

As a concluding remark we now discuss our result 
in a more general context. Early QMC simulationP^ 
on clean 2DEGs at R = 0 T showed that the en¬ 
ergy of the spin and valley polarised system becomes 
lower than that of the two component system in 
the region Vg > 20 before Wigner crystallisation at 
Tg ~ 34. A ferromagnetic instability had therefore 
been expected in the valley polarised system. In 
contrast, and in qualitative contradiction with re¬ 
sults from Hartree-Fock calculations, no such insta¬ 
bility was observed in the QMC simulations for the 
valley and spin degenerate system which was found 
to be the stable phase all the way up to Wigner 
crystallisation at Vg 42 (See Fig. 1 in Ref. HU) . 
Our experiments and quantum Monte Carlo simu¬ 
lations confirm and demonstrate that this scenario 
remains valid in real disordered systems, rather than 
it being only a special case of hypothetical disorder- 
free systems. Results from simulations displayed 
in Fig. show that Rp(pv = I) is approaching 0 
at a finite electron density indicating a ferromag¬ 
netic instability in the valley polarized system. The 
instability occurs at a lower Vg in our disordered 



Figure 6. Density dependence of Bp. Theoretical 
dependence of Bp in a single particle pictnre and in pres¬ 
ence of interactions for various values of disorder (quan¬ 
tum Monte-Carlo simulation). 


system compared to the clean onJ^. In contrast, 
no sign of such instability is seen in the curve for 
Bp{pv = 0) down to the lowest electron densities 
we have studied. Experimentally, the mobility of 
our sample obviously does not allow us to reach the 
very low densities necessary for the observation of 
spontaneous spin polarisation or Wigner crystallisa¬ 
tion. Yet, the crossing of the curves Rp(pv = 0) and 
Bp(jpv = 1) provides strong experimental evidence 
of the higher stability of the spin-valley degenerate 
system because they demonstrate that interaction- 
induced spin alignment is much less efficient in the 
valley unpolarised system than in the valley po¬ 
larised system. The strong enhancement of spin 
susceptibility with valley polarisation in AIAJ^^ also 
supports this interpretation. The excellent agree¬ 
ment between the theory and our experiments sug¬ 
gests that the result can be extrapolated to cleaner 
spin-valley degenerate systems with greater inter¬ 
action. In those systems, in the absence of ferro¬ 
magnetic instability, we anticipate the presence of a 
strongly correlated electron liquid. This may result 
in a rich physics which might soon be accessed ex- 
ploitin g rece nt developments in high mobility silicon 
systemP^. 
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METHODS 

Samples and control of valley polarisation 

The samples consists of a Si02/Si(100)/Si02 
quantum well of nominally lOnm thick silicon with 
front- and back-gate oxide thicknesses of 75 nm and 
380 nm respectively. The fabrication proceedure of 
this sample is described in Ref. [22l A degenerately 
phosphorus-doped polysilicon layer was used as front 
gate while the substrate was used as a back gate. 

It is well known that valley degeneracy can be 
lifted at the (001) Si/Si02 interface in Si MOSFETs 
and that the valley splitting can be increased by in¬ 
creasing the out-of-plane electric held, which in tra¬ 
ditional MOSFETs can be controlled by changing 
the substrate biaJi^. The magnitude of this valley 
splitting is found also to depend on the way in which 
the Si-Si 02 interface is prepared, and the use of a 
buried-oxide interface using SIMOX Reparation by 
IMplantation of Oxygen) technology^ allows us to 
enhance the valley splitting up to tens of meV^. 
The coupling responsible for the bare single-particle 
splitting is indu ced pr edominantly by the large inter¬ 
face electric helcP^Hll but there remain uncertainties 
as to the exact microscopic details that give rise to 
the particularly large values in SIMOX buried-oxide 
interface^^^l*^. 

Experimentally, the valleys splitting is determined 
by htting Shubnikov de Haas oscillations with an 
empirical expression for valley splittin^^ 

Av = a6n (2) 

where, Sn is an empirical measure of the out-of- 
plane electrostatic potential asymmetry controlled 
by front and back gates: 

5n = ub — n-p (3) 

where np and np are electron densities contributed 
by respective gates. Both np and ub can take pos¬ 
itive or negative values where negative values rep¬ 
resent a density reduction due to a depleting bias 
from the corresponding gate so that the total elec¬ 
tron density is given by 

n = Ub + Up. (4) 

The numerical factor a which is of the order of 0.5 
meV/10^^m“^ determines how much the valley split¬ 
ting changes with Sn when Sn is positive. That is, 
when the quantum well is biased in such a way that 
the electrons are pulled towards the back (SIMOX 


buried-oxide) interface. On the other hand, when 
Sn is negative, the electrons are pushed against the 
front interface, which is formed by standard ther¬ 
mal oxidation, where we hnd the valley splitting to 
be negligibly small. Thus, by pressing the electrons 
against the buried-oxide interface (positive Sn) we 
can increase the valley splitting continuously, and 
independently control the electron density n. The 
out-of-plane potential necessarily affects the disor¬ 
der, however, but the effects of this can be indepen¬ 
dently examined by applying a negative Sn for which 
there is no valley splittingAl_ 

For fitting the Shubnikov de Haas oscillations, we 
fix the perpendicular magnetic field and compare the 
valley splitting Av against the cyclotron energy HoJc, 
or more accurately, we map the number of occu¬ 
pied Landau levels of the two valleys as function 
of {Sn,n). Straightforwardly applying the single¬ 
particle model only yields a value for amp but not a, 
in the same manner as coincidence experiments un¬ 
der tilted field only provide values for gm^ and not 
g. The bare effective mass only provides a crude con¬ 
version of the valley splitting to an energy scal^^H^ 
and would represent a good measure of the valley 
splitting in the absence of interactions that alter the 
effective density of states. 

Valley polarisation is then estimated at full spin 
polarisation from the equation: 

AyDo 


The determination of valley polarisation does not, 
therefore, rely on separating a and ms since in 
Eqs. Av never appears on its own but always as 
a product with Dq. It follows that this equation re¬ 
mains a valid method of determining the valley po¬ 
larisation even in the presence of strong interactions. 
See Supplementary Discussion for more details. 


Electrical Measurements 

The samples were cooled in a Variable Temper¬ 
ature Insert with a base temperature of 1.6 K in¬ 
serted into a 30 T resistive magnet. A standard 
four-terminal lock-in technique was used to measure 
the resistivity pxx and pxy of a sample with a Hall- 
bar geometry. The current was kept below 5 nA to 
avoid electron heating. The samples were aligned 
parallel to the applied magnetic field with an in-situ 
rotator, eliminating the Hall resistance. 





Determination of the mobility 

The comparison of the theory to the experimen¬ 
tal data requires the determination of the mobil¬ 
ity of the sample in the Drude regime (conductiv¬ 
ity >> e^/h, see the next section). In practice one 
should measure the resistance at high density to ob¬ 
tain this quantity. However, this is not straightfor¬ 
ward in our case because for densities larger than 
n = 5 X 10^® m“^, electrons start to experience scat¬ 
tering from localized states in the upper spatial sub¬ 
band which become populatecPSl. Fortunately, the 
mobility at high density in the absence of the influ¬ 
ence of the upper spatial subband can still be de¬ 
termined in our sample. To do so, we exploit the 
fact that pressing the electron gas to the front inter¬ 
face at negative 5n does not increase valley splitting 
but narrows the effective width of the out-of-plane 
electronic wavefunction. This pushes the upper spa¬ 
tial subband to higher energies and suppresses its 
influence. This is illustrated by the solid white line 
in Fig. which demarcates the onset of occupation 
of the upper subband. As the out-of-plane elec¬ 
tric bias {5n) is increased, the density n at which 
it starts to fill increases, reflecting the increasing 
confinement energjl^. The suppression of scatter¬ 
ing by localized states in the upper spatial subband 
is evidenced by the reduction of resistance with 5n 
for densities about 10^® m“^. For densities relevant 
to the present study, scattering by localized states 
in the upper subband is absent as shown by the 
blue constant resistance contour in Fig. which is 
parallel to the constant density line for all negative 
—2 X < dn < 0. For even larger electric 

fields 6n < —2 x the electron g^ experi¬ 

ences the roughness of the front interface^ and the 
resistance starts to increase so that the iso-resistance 
lines deviate from constant density lines. We con¬ 
clude that the mobility at high density and in ab¬ 
sence of the influence of the upper spatial subband 
can be estimated in the regime —2 x < 

i5n < — 1 X 10^®m“^. Therefore, we estimated ^ for 
n = 10^® m“^ and Sn = —1.5 x marked 

by a star in the Fig. where p = 760H leading to 
p, = 8000 cm^V“^s“^. 


Quantum Monte Carlo simulations 

The model we consider is a generalization of the 
Anderson model to the many-body problem (see 
Ref. HI] and ED for details). The system is made of 
N spin up/down electrons with Coulomb interaction 


on a disordered lattice of x Ly sites. The elec¬ 
trons either populate a single valley (pv = 1) or are 
equally split up into two degenerate valleys (pv = 0). 
Their spin conhguration sets the value of the spin 
polarisation ps. Formally, the spin and valley de¬ 
grees of freedom are treated strictly in the same way 
as an internal electronic degree of freedom. In the 
continuum limit v = N/L^Ly ^ 1 (where lattice ef¬ 
fects are negligible) and in the thermodynamics limit 
1, the physics of a given (ps,pv) configuration 
is entirely controlled by two dimensionless parame¬ 
ters Ts = rrihe^/{ATreti^y/im) (mb = 0.19me effective 
mass, e electron charge, e = 7.7eo dielectric con¬ 
stant) and l/kpl {kp Fermi momentum, I mean free 
path, both taken for the spin- and valley-degenerate 
system) which characterize respectively the interac¬ 
tion strength and the disorder strength in the sys¬ 
tem. Experimentally, it is difficult to estimate kpl 
for low density systems because conductivity is no 
longer a good estimate of disorder. We overcome 
this issued by observing that for white noise dis¬ 
order of amplitude W (as assumed in our model), 
kpl oc n. Therefore, p = rsy/kpl does not depend 
on n and depends only on H®1. One can estimate 
rj in the high density regime where electronic inter¬ 
actions are negligible. In that regime, the conduc¬ 
tance g of the system is p = {2e^/h)kpl which gives 
p = ^/pe®/^m/(47re/i®/^) where p = gl{en) is the 
mobility of the sample. Thus, the two input param¬ 
eters Tg and kpl of our model can be estimated from 
the mobility of the sample measured at high den¬ 
sity and from the experimental values of electronic 
densities. 

We use the Green’s Function Monte Carlo 
methocP^ in the fixed-node approximation to 
compute the energy per particle E(ps) of the 
ground state of our model at zero temperatur^Sl 
The polarisation energy Ep is deduced from 
Ep = E{ps = 1) — E{ps = 0) and averaged over 50 
to 200 samples depending on the disorder strength. 
To extrapolate data at the thermodynamic and 
continuum limit, finite N- and z^-effects are carefully 
investigated. We thus observed large but controlled 
lattice effects without interaction that disappear as 
interaction is switched on (rg > 0.5). Small finite 
size effects in N are also present but they rapidly 
fade with the disorder amplitude. The resulting 
extrapolated data, as well as their fits given below, 
are finally obtained with a precision of the order 
of ±0.02F;0 (for kpo > 1-5) and iO.OdE® (for 
0.3 < kplo < 1.5) at pv = 0, and roughly twice 
larger at Pv = 1- 

Without interaction (rg = 0), our data are in per- 
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feet agreement with the second-order perturbative 
formula Ep/E^ = 1/2 [ 1 ] -|- log 2/(7rfc^Zo) for Pv = 0 
[Pv = l]j at least for weak disorder ^ 0-4) ■ 

In the presence of (even small) interaction, first the 
effect of disorder is reversed making easier the spin 
polarisation of the system and second, the 1/kPplQ- 
correction to Ep is no longer valid (except for tiny 
disorder). We find that our Ep data are very well 
described by the following formula, 

Ep{k%l^, r,) = E;\t,) + + ^5 (6) 

\/kplo 

where the polarisation energy of the clean system 
Ep and the parameter /3 are both fitted with Fade 
approximates. 


and interaction strength (0 < rg < 10 , > !)■ In 

particular, Eq. is satisfied in the disorder and in¬ 
teraction regime explored in the present experiment. 
Then, when an in-plane magnetic field B is applied, 
a Zeeman term —p/iB^Ps/S has to be added to the 
right hand side of Eq. Minimizing the energy 
E with respect to Ps gives the spin polarisation of 
the system at zero temperature, /[AEp), from 

which we get Bp = AEp/{gg,p,)- 
Figure presents the result of the numerical cal¬ 
culations of the magnetic field of full spin polarisa¬ 
tion for various values of disorder. This figure shows 
that even in the presence of weak disorder (p = 10 ® 
cm^V~^s“^), the magnetic field of full spin polar¬ 
isation is expected to be much lower than in the 
non-interacting picture. 


Tpclf„ \ _ ^0 + ^1 fs J^O 

- A, + A,r. + A,rl <"> 


/?(f’s) = 


±2 -h ^3 '/’s -r ^4 ^ 

Bq + Bi rl 
B 2 -b ^3 Tg -b B 4 r2 ’ 


( 8 ) 


the fitting parameters Ai and Bi being given in Ta¬ 
ble 1^ We note that Eq. for Ep is in very good 
agreement with previous Quantum Monte Carlo 
calculations performed for the valley-polarized sys- 
tenP^ and the valley-degenerate systenP^. Equa¬ 
tion]^- and in particular the fact that /3’s sign flips 
at rg « 0.3 (for pv = 0) and Tg « 1.2 (for pv = 1) 
- mainly depicts the opposite effects of disorder at 
very weak and stronger interactions. The last pa¬ 
rameter adjusts the origin of the linear disorder 
correction in 1 / \/kplQ, to roughly take into account 
the actual quadratic correction in 1 /\/kplQ at weak 
disorder. At extremely weak disorder, l/{kpP) < 
0.04, A 5 shall be taken equal t(^^and Ep evaluated 


by Ep. We point out that Eqs. ^ 


are no more than 
one simple way to report our dafa, valid (at least) 
for 0.25 [0.5] < Tg < 10 (at Pv = 0 [pv = Ij) and 
k^plo > 0.3 as long as the output Ep is positive. 

Deducing the polarisation magnetic field Bp from 
Ep is straightforward, once noticing that the energy 
E of the ground state is quadratic in ps, 


E(pg) = E{0) + Eppl 


(9) 


This statement is obvious in the absence of disorder 
and interaction where we have Ep = Ep/2 [Ep] for 
Pv = 0 [pv = 1]. Numerically, it turns out to remain 
valid with good precision for intermediate disorder 


i 

0 

1 

2 

3 

4 

5 

o' 

II 

> 

27.93 

9.83 

56.5 

46 

1.77 

0.019 

Aiijp^ = 1 ) 

23.4 

-0.5 

24.4 

7.75 

0.27 

0 

o' 

II 

> 

2.70 

-24.81 

42.4 

-11.5 

246.9 


Bi(Pv = 1 ) 

54.28 

-39.30 

348 

170.5 

267 



TABLE I. Coefficients for the estimation of polari¬ 
sation energies. Ai and Bi parameters of Eqs. [^and|^ 
for the valley-degenerate (pv = 0) and valley-polarized 
(pv = 1) system. 

The code to perform this simulations has been 
parallelised and ported on CEA Computing Center 
for Research and Technology (CCRT) massive par¬ 
allel clusters. About 100000 CPU hours have been 
required for the present study. The code is available 
upon request. 


TEM images 

The High Angle Annular Dark Field scan¬ 
ning TEM images were measured using a probe- 
aberration corrected FEI Titan microscope operated 
at 200 kV. A 100-nm-thick specimen was prepared 
by focused ion beam milling at 5 kV to reduce the 
surface damage. The HAADF STEM images are 
sensitive to Z-contrast and the vertical bright dumb¬ 
bell structures are typical of aberration corrected 
images of silicon samples oriented in the < 110 > 
direction and show silicon atoms separated by 1.36 
Angstroms in projection. 
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